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1. Introduction 


At the request of and with support from the Naval Research Laboratory*, a 
comprehensive study was conducted to determine the space radiation environment 
of a series of high inclination geocentric orbits for a mission duration of 
five years. | 


‘SFollowing the precedent established with—previous—seudfes, the external 
(surface incident) charged particle radiation, predicted for the satellite was 


determined by orbital flux integration for six independent trajectories, (see 
seetion 3}> The latest standard models of the environment were used in the 
calculations <tsee-pectron 5). 


‘>Magnetic field definitions for the six nominal circular trajectories were 
obtained from a current field model, (see-section 3). 


‘Spatial and temporal variations or conditions affecting the static environment 
models were considered and accounted for, wherever possible. 


‘Finally, limited shielding and dose evaluations were performed for simple 
infinite slab and spherical geometries. 


Results, given in graphical and tabular form, are analyzed, explained, and 
discussed. Conclusions are presented and commented on. / 


2. SPECIFICATION OF ORBITS 


The analysis was based on nominal circular orbits with inclination of 60 
degrees and altitudes of 1667, 2593, 3889, 5186, 6389, and 10371 kilometers. 


3. 


Six separate flight path emphemerides were generated for the specified orbits 
with the GEODYN-BLCONV System! for trajectories of 24-hour duration defined at 
i-minute intervals. The length of simulated orbit time and the integration 
stepsize were especially selected so as to provide sufficient point density to 
insure an adequate sampling of the ambient radiation environment when flying 
the trajectories through the models. The trajectories were subsequently 
converted from geodetic polar to magnetic B-L coordinates with McIlwain's 
INVAR program of 19652 and the field routine ALLMAG,? which now utilizes the 
BARRACLOUGH 1975 field model.4 The field computations were extrapolated to 
the tentative mission epoch of 1989.5 with linear time terms representing 
secular variations of the field. 


4. FLIGHT PATH EXPOSURE TO TRAPPING DOMAINS 


The specified nominal flight-path configurations display a significant 
characteristic of high inclination orbits in magnetic L-space: they traverse 
almost the entire terrestrial radiation belt twice during each revolution, 
moving back and forth through regions of low L values (the inner zone: 1.0 < 
L < 2.8), regions of high L values (the outer zone: 2.8 < L < 12), and 


*This work funded by the Naval Electronics Systems Command under the NRL 
Nuclear Survivability/Vulnerability Program. 
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regions outside the trapping domain (external). Occasionally, some 
revolutions will also enter regions of space where no particle trapping can 
occur because of atmospheric cut-off conditions; that is, trajectory segments 
may have a combination of magnetic B and L values that place them outside the 
atmospheric cut-off limits of the models. 


These excursions and the “external” visitations afford the satellite an amount 
of flux-free time, which may be of substantial duration (see section 11, C). 


5. TRAPPED PARTICLE ENVIRONMENT MODELS 


The fluxes in this study were obtained from current NSSDC models: the solar 

maximum AE6 for the inner zone electrons5, the new interim model AEI7 for the ° 
outer zone electrons®, and the solar maximum version of the new AP8 model’ for 
energetic trapped protons. It should be noted that the interim AEI7 does not 
reflect solar cycle variations in its present state. However, this model was 
issued in two versions, the AEI7-HI and the AEI7-LO, in order to account for 
differences in the data sets used in their construction. The LO version was 
used in this effort. All models describe an average static environment at a 
given epoch. 


6. ORBITAL FLUX INTEGRATIONS 


Orbital flux integrations were performed with the UNIFLUX® and the SOFIP? 
systems. UNIFLUX provides L-band distributions and exposure times with B-L 
bin breakdown, while SOFIP provides the dose and shield data. 


7. GEOMAGNETIC SHIELDING AND SOLAR FLARE PROTONS 


Low altitude high inclination orbits experience a significant amount of 
geomagnetic shielding from cosmic rays of solar or galactic origin in the 
energy range E > 10 MeV. Therefore, it may be assumed that the spacecraft 
will only intermittently be exposed to the unattenuated interplanetary solar 
flare proton intensities of all energies above 10 MeV. To a first 
approximation, the fluxes may also be considered omnidirectional and 
isotropic, probably to within 10-15%. 


Usually, geomagnetic shielding effects on geocentric missions are being 
evaluated with simple rigidity considerations because of substantial diurnal 
variations in the cutoff latitude associated with geomagnetic tail effects 
(2-4 degrees) and storm-induced changes (> 4 degrees). The simple analysis 
used here assumed that energetic solar protons of all energies above 10 MeV 
have free access to all magnetospheric regions external to a dipole shell of 
L*5 earth radii, which is equivalent to a cut-off latitude of about 63 
degrees. as 


Predictions of solar flare proton fluxes at 1 AU are obtained as a function of 
mission duration 1 and confidence level 9* on the basis of a probabilistic 
analysis!0 using a modified type of Poisson statistics by a computerized model 
SOLPRO'! that includes the distinction between "ordinary" (OR) and 


——— 
*Q denotes the degree of confidence one wishes to assign to the results, namely 


that for the specified mission duration the calculated fluences are the 
smallest values which will not be exceeded by actually encountered intensities. 


2 


"anomalously large" (AL) events and the probability of occurrence of the 
latter. Both AL=- and OR- event fluences are non-linear functions of Q and T. 
For these predictions, only high quality comprehensive satellite measurements 
(not ground observations) are being used, covering almost the entire 20th 
solar cycle. There have been indications that descriptions of the solar flare 
environment in interplanetary space (at 1 AU), derived from interpretations 
and extrapolations of ground based measurements, have not been very accurate. 


It should be noted that the statistics cannot predict when an AL event will 
occur; only the probability that one will occur in a given length of time. 

And it must be remembered that a single AL event will impart its total fluence 
within two to four days. 


This implies that for unmanned satellites with mission durations of T >? 1 
year, OR-event fluences are not significant because probabilistic theory 
predicts the possible occurrence of at least one AL event, even for the lowest 
allowable confidence level (Q=80%). 


8. FLUX DATA: TYPE, QUALITY, AND VARIATIONS 


The trapped particle flux data available from the models represent 
omnidirectional, integral intensities that one would expect to obtain as 
average values over periods in excess of six months. But over most regions of 
magnetospheric space (L>2 earth radii), short term excursions can vary from 
these values by factors of 102 to 103, depending on the particle energies and 
on the type and intensity of the causative event. These variations do affect 
the investigated missions because their trajectories enter regions of space 
where L is greater than 2 earth radii. Also, trapped particle populations 
experience changes due to: (a) local time (LT) dependence, and (b) solar cycle 
dependence. Both are of some consequence to these missions. The former is 
significant for spacecraft that sample regions of L > 5 earth radii, which are 
visited by the specified missions. To compensate for these variations, the 
model provides LT-averaged values, which should yield an adequate approximation 
for missions of long duration (t > 1 year). The solar cycle variations have 
been taken into account by selecting the appropriate models for each period, 
where available. 


Generally, solar cycle variations have opposite effects on each particle 
species: 


Solar Min Solar Max 
Electrons lower higher 
Protons higher lower 


The solar cycle changes, as derived from a comparison of the corresponding 
models, are functions of energy E and magnetic parameter L. For the inner 
zone electrons, they may range from a factor of 1 to a factor of 5. 


Protons are only affected in the vicinity of the atmospheric cutoff regions. 
No changes of consequence have been observed in the heart of the proton 
trapping domain. Proton changes have about the same range as those of the 
electrons. 


It is necessary to emphasize that the calculations, although based on the best 
data available for the past epochs, can only serve as approximations for the 
future. 


It also should be noted that a basic uncertainty factor of 2 is defined for 
the flux values of the AP8 and the AE5 models, while the AE6 is characterized 
by an average uncertainty factor of 5. No uncertainty factor has yet been 
defined for the interim AEI7. 


' 9. DOSE AND SHIELDING EVALUATION 


Doses were calculated from the total orbit integrated, surface incident, 
omnidirectional, integral fluences by existing shielding codes!2, as functions i 
of various aluminum shield thicknesses and geometries. 


A simple procedure was followed, not involving solid angle sectoring or 

three-dimensional ray tracing considerations. Instead, a simple 

two-dimensional geometry with a cosine law for the incident spectra, and a 
three-dimensional spherical geometry were considered. (See comment in section 

11D-ITI) 


Bremsstrahlung calculations were performed with the same codes. 
10. RESULTS: PRESENTATION DESCRIPTION 


This section describes the form and format in which the results, derived from 
the Orbital Flux Integration (OFI) process, are presented for practical use. 

Except where otherwise specified, all particle data in this report relate to 

integral, omnidirectional fluxes or fluences. 


A. Tabular Presentations 


The outcome of all calculations is summarized in Tables 1 to 60. The tables 
are arranged in six sets, where every set pertains to one specific type of 
data. The first two sets have two similar members for every trajectory 
considered in the study: one for trapped protons and one for electrons, in 
that order. The next three sets contain only one member for each trajectory. 
The sixth set contains tnree similar members for each trajectory. A more 
detailed description of the tables is provided in the following paragraphs. 


I. L-band Tabulations: Tables 1-12 
Tabulation of total orbit-integrated fluence distributions by L-bands for 


selected energy thresholds, in units of particles per square centimeter, 
normalized to 5 MeV and .5 MeV for protons and electrons, respectively. « 
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: The tables contain 48 L-bands of equal size covering the range from L=1.0 to 
- L > 10.4 earth radii in constant increments of .2 earth radii. 
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II. Spectral Profiles: Tables 13-24 


a 
? Tabulation of average orbit-integrated spectral distributions. Composite 
spectra are given in units of: fluxes per square centimeter per second, 


| fluxes per square centimeter per day, and total fluences per specified mission 
i ft duration (5 years). For the electrons, the latter are also given in terms of 


: 
in 


inner and outer zone contributions. Functionally derived differential fluxes 
are listed in the last columns for both species of particles. 


Total orbit-integrated spectra in percent, for energy intervals AE 
corresponding to the energy levels of the L-band tables, are also given in 
terms of average instantaneous and daily intensities. 


An exposure index (for the normalization energies used in the L-band tables) 
is listed for nine successive intensity ranges varying by one order of 
magnitude, in terms of processed exposure duration (in hours) and total number 
of particles accumulated while in that intensity range for the indicated 
number of hours. 


III. Peaks and Totals Per Orbit: Tables 25-30 


These tables contain the absolute instantaneous peak fluxes and the total 
fluences accumulated during each successive revolution, as obtained from 
the nominal trajectories for the investigated flight duration (24 hours of 
mission time). 


Specifically, there are nine columns on these tables. Column 1 is an orbit 
counting device, based on: 
a) the orbit period when the trajectory is circular and lies in the 
equatorial plane; 
b) the physical perigee in all elliptical flight-path cases; and, 
c) the equatorial crossing for circular inclined trajectories. 


Column 2 gives the peak flux. Columns 3, 4 and 5 indicate the spacecraft 
position in geocentric coordinates at which the predicted peak flux was 
encountered. Columns 6, 7 and 8 determine respectively the relative orbit 
time and the magnetic B-L coordinates for this event. For the purpose of 
orbital radiation studies, all simulated trajectories start at t, = hours. 
Finally, the last column indicates the total predicted flux to he encountered 
during that particular orbit. It is advisable to disregard the last line on 
this table because many times that orbit is incomplete and the fluxes or 
positions shown do not correspond to true peaks. 


Iv. Time-Accounting and Exposure-Analysis: Tables 31-36 


The “EXPOSURE-ANALYSIS" summary indicates what percent of its total lifetime T 
the satellite spends in “flux-free" regions of space, what percent of its 
total lifetime it spends in high intensity proton and electron domains, and 
while so exposed, what percent of its total flux it accumulates. 


In the context of this study, the term "flux-free" applies to all regions of 
space where trapped particle fluxes are less than one proton or electron per 
square centimeter per second, having energies E > 5 MeV, and E > .5 MeV, 
respectively. By definition, this includes all regions external to the Van 
Allen radiation belts. 


The concept of “trapped particle fluxes" is meant to include stably trapped, 
pseudo trapped, and transient fluxes, as long as they are part of or contained 
in the environment models used and, in the case of transients or pseudos, 
their sources are considered powerful enough to supply them continuously in 
substantial numbers. 


Similarly, as "high intensity" are defined those regions of space where the 
instantaneous, integral, omnidirectional, trapped-particle flux is greater 
than 103 protons with energies E > 5. MeV, and greater than 10° electrons with 
energies E > 5. Mev. 


The values given in these tables are statistical averages, obtained over 
extended intervals of mission time. However, they may vary significantly from 
one orbit to the next, when individual revolutions are considered. 


The “TIME-ACCOUNT" breakdown shows what percent of its total time the 
satellite spends in the “inner zone" (1.0 < L < 2.8) and in the “outer zone" 
(2.8 < L <¢ 11.0) electron trapping domains, and also the percent of time spent 
in regions external to the latter (L > 11.0). 


It should be noted that the confinement of the outer zone within the boundary 
of the L=11.0 earth radii volume is arbitrary and has no physical meaning. It 
is intended only as a simplification to facilitate the calculations. The 
region considered “external” in this study (L > 11.0) is still partially a 
domain of the outer zone, at least as far out as L=12.9 earth radii, according 
to the current environment models. 


A last item on this table: the inner zone time is further subdivided into two 
parts: the percentage of time spent outside (L < 1.1) and inside (1.1 < L < 
2.8) the trapping domain. 


Vv. Solar Proton Fluences and Exposure Factor: Tables 37-42 


For the specified mission duration T (printed in the sub-title), and dipole 
cut-off shell (L=5 earth radii, shown in the header), this table lists the 
solar proton fluence-spectra (in units of particles per square centimeter) at 
five discrete confidence levels 9 (given at the top of each column). 


The exposure factors (in percent of total mission duration) obtained from the 
geomagnetic shielding analysis are also listed for four dipole cut-off shell 
values (in earth radii). 


VI. Total Dose and Components: Tables 43-60 


These tables list doses in units of radsay as a function of aluminum shield 
thickness, given in three ways: range s in grams per square centimeter, depth 
t in millimeters, and depth t in mils. 


Electron, bremsstrahlung, and proton contributions to the overall sumtotal 

dose are given separately. Electron and proton doses are further broken down 

into their respective constituents; namely, inner zone and outer zone for the é 
former, trapped and solar flare for the latter. 


The specific mission duration for which the doses have been calculated is 
indicated in the table headline. 


Caution: the AL-event solar flare protons are not contributed gradually over 
the investigated mission duration (t = 5 years) but are imparted in toto in a 
relatively short burst, that is, within approximately 2-4 days per AL event. 


B. Graphical Presentation 


Some of the tabulated data are also plotted in Figures 1 to 36, and 49 to 102 
with additional Figures 37 to 48 containing plots of flight path data. 
Positional flux and dose data are plotted in Figures 103-162. As with the 
tables, the computer plots are arranged in eight sets, where again each set 
pertains to one specific type of data. The first three sets have two similar 
members for each trajectory investigated: one for each particle species. The 
next two sets (fourth and fifth) contain one member for each trajectory 
considered. The sixth set contains nine similar members for each trajectory, 
providing three graphs (for respective depth ranges) for each of three 
geometries. The seventh set contains two similar members for each trajectory: 
one for each particle species. Finally, the last set has eight members for 
each trajectory. 


I. Time and Flux Histograms: Figures 1-12 


These plots show two curves superimposed on the same graph; rely, one each 
for the variables "time" and “flux". Both are given on a: ‘-log scale as 
functions of the parameter L (earth radii), within the ran 1<L< 10, and 
for constant L-bands of .1 earth radius width. The plots c p‘ 

a) by a plain curve, the characteristic trajectory inte -v1ties as 
obtained from the orbital integration process in terms of averaged 
integral particle fluxes above a given energy. 

b) by a contour marked with suymbols, the percent total lifetime (%T) 
spent in each L interval. 


The logarithmic ordinate relates to the time-flux variables. The printed 
numbers are powers of 10 and pertain to the fluxes; the scale values for the 
time curves are given in the upper part of the ordinate label: from 1073 to 
102 percent of T, the type of particles, their integral energy, and the units, 
are all given in the lower part of the label. The label on the top of the 
graph identifies the trajectory. 


II. Spectral Profiles: Figures 13-24 


A graphical presentation of the final composite spectral distribution, 
obtained from the orbital integration process. The plots are semi-log graphs, 
where the abscissa is a linear energy scale for integral particle energies E, 
in MeV, and the ordinate is a logarithmic scale for the fluxes, given in 
daily averages for energies greater than E; the printed scale values are 
powers of 10. 


III. Peaks Per Orbit: Figures 25-36 


Here the absolute peak intensities, encountered per period (1 period = 1 
revolution = 1 orbit), are plotted for the duration of the flight-time 
processed in the analysis. The logarithmic ordinate, with scale values in 
powers of 10, relates to instantaneous particle fluxes of the environment at 
the indicated energy thresholds, while the abscissa is a linear orbit 
enumeration. 


Iv. Trajectory World Map Projections: Figures 37-42 


These graphs depict the surface trace of the geocentrically projected 
subsatellite positions. The trajectories are plotted for several revolutions 
on a qlobal map produced by a Miller Cylindrical Projection method. The 
contours of the continents have heen omitted for clarity. The positions of 
equatorial crossing, of physical perigee, or of period commencement are 
indicated by numhers identifying the orhits shown in the graphs. For all 
trajectories, the distance between successive sequential numbers is a measure 
of the orbit precession. 


Vv. Flight Path Tracing in B-L Space: Figures 43-48 


Plots showing trajectory traces in B-L space on a semi-log scale. Several 
orbits are depicted, each identified by its sequential number. The magnetic 
equator is entered on all plots. The logarithmic ordinate relates to the 
field strength B in gauss; the printed values are exponents of 10. L is given 
in earth radii on the linear abscissa. 


VI. Dose-Depth Curves by Geometry: Figures 49-102 


Plots of final depth-dose values for the indicated mission duration. 
Normally, these plots show composite curves for bremsstrahlung, combined 
electrons (inner and outer zone), combined protons (trapped and solar flare), 
and sumtotal of all contributions. In the present case, the respective 
contours consist of inner and outer zone electrons separately and of trapped 
and solar flare protons, separately, of composite bremsstrahlung, and the 
sumtotal. 


For ease of use and in order to nrovide a greater .esotution at the more 
sensitive range of depths (namely the thinner shields) three plots have been 
generated per processed trajectory, for shield-ranges and subdivisions 
increasing by one order of magnitude. 


The logarithmic ordinate, with scale values in powers of 10, relates to 
aluminum dose in units of rads. The linear abscissa is the shield thickness, 
given in three different units: range s in grams per square centimeter, depth 
t in millimeters, depth t in mils. 


VII. Positional Flux Plots: Figures 103-114 


Plots of instantaneous omnidirectional trapped particle fluxes (electrons and 
protons) at (up to 10) specified threshold energy levels (>MeV), for a 
selected orbit (usually worst case revolution through heart of SAA). 


The logarithmic ordinate, with scale values in powers of 10, relates to the 
number of particles per square centimeter per second. The linear abscissa is 
the relative time, in minutes or hours , from the beginning of the selected 
orbital pass. 
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VIII. Positional Dose Plots: Figures 115-162 


Plots of instantaneous omnidirectional trapped particle dose values at (up to 
10) specified shield thicknesses (omnidirectional isotropic incidence, 
cosine-theta distribution) for a selected orbit (usually worst case revolution 
through heart of SAA). Separate plots are generated (if present) for: 

electron dose (including bremsstrahlung), proton dose, and total dose (no solar 
proton contributions are included) for dose at transmission surface of aluminum 
slab shields, dose in semi-infinite aluminum medium, and dose at center of 
aluminum spheres. 


The logarithmic ordinate, with scale values in powers of 10, relates to the 
respective dose in units of rads-aluminum. The linear abscissa is the relative 
time, in minutes or hours, from the heginning of the selected orbital pass. 


11. RESULTS: ANALYSIS AND DISCUSSION 


In this section, some of the presented tabular or graphical study~-results are 
discussed, with occasional comments as to their use, limits, and 
applications. 


A. Spectral Profiles 


Characteristic features of the near earth radiation environment are strong 
altitude and inclination dependencies. However, at high inclination values 
(30° < i < 90°) small changes in inclination will not produce changes in flux 
levels and spectral distributions as significant as those produced by small 
changes in altitude. The greatest inclination dependent variations occur in 
the range 0° < i ¢ 30°. 


I. Protons: The protons exhibit relatively hard, almost uniform 
spectra in the investigated inclination-altitude regime. 


It should be noted that a characteristic softening of the high energy tail of 
the spectrum appears at the 6389 km altitude orbit, near E > 400 MeV, as a 
consequence of the limit of the volume of space occupied by these particles. 
The extent of the proton trapping domain along the maqnetic equator is 
inversly related to their energy. This is strikingly demonstrated by the 
spectrum of the 10371 km altitude orbit, which falls off sharply at about 

E > 100 MeV and which indicates that according to the standard model there are 
no protons with energies E > 150 MeV contained in the region of space sampled 
by this trajectory. 


II. Electrons: The electrons show complex variable spectra. 
Inner zone and outer zone average, orbit-integrated, composite intensities 
rise non-uniformly with altitude, particularly at energies above 3.75 MeV with 
differences reaching up to several orders of magnitude at E > 6 MeV. Spectra 
extend to higher energies as height increases. 


The inner zone spectra fall rapidly off to zero flux in the energy range from 
4 to 5 MeV and they are therefore more benign then their harder outer-zone 
counterparts, which extend to energies of about 7 MeV. 
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At low altitudes (1667 abd 2593 km), the inner zone contributions prevail up 
to about E > 2.75 MeV, but for all energies E > 2 MeV the inner zone spectra 
are always softer than the outer zone spectra (in agreement with the modeis). 


However, at high altitudes (e.g. the 10371 km orhit), the inner zone 
contributions become insignificant even at low energies (.1 < E(MeV) < .7). 


B. Peaks Per Orbit 


The absolute peaks per revolution have been obtained for standard processing 
energies: E> 5. MeV for protons and E > .5 MeV for electrons. Other energy 
selections produce different peak curves in an inverse relationship: lower 
energies yield higher and more expanded contours, and vice versa. 


Peak contours of inclined circular trajectories display amplitude variations 
and sometimes discontinuities (flux-free time) that follow periodic patterns 
based on the daily cycle of revolutions. For fixed energies, amplitudes and 
discontinuities are function of: (a) inclination i, and (b) altitude h. 


Variations in either i or h may produce significant changes in the amplitude 
of the peak curves and in the duration of the discontinuities: up to several 
orders of magnitude for the former, and completely eliminating the latter. 


For the investigated trajectories at the given, fixed, inclination (i = 60°), 
the following observations can be made: 

a) protons: at the processed energy of E > 5 MeV only small variations 
are obtained for the low altitude orbits at 1667 and 2593 km; the difference 
between maximum and minimum peak intensity predicted is about a factor of 6 
ar 13, respectively. These differences disappear at higher altitudes, while 
at the sade time the peak flux levels rise by about two orders of magnitude, 
from ~7x104 particles per square centimeter at 1667 km to ~7x10 particles per 
square centimeter at 5186 km. 


b) electrons: at the processed energy of E > .5 MeV the electron peaks 
appear very similar to those of the protons, as described in the previous 
paragraph. But the electron peak fluxes start at higher levels and reach an 
upper limit at the same altitude as the protons only with smaller differences: 
from ~3x10© particles per square centimeter at 1667 km to ~1x107 particles per 
square centimeter at 5186 km. 


Cc. Flux-Free Time 


Some comments on this topic have been provided in the previous section and in 
section 10/IV. Here a more detailed discussion will be given. 


Flux-free time (FFT) intervals are an important feature of certain orbital 
configurations. They may occur over short orbit segments (partial FFT per 
period) or over the entire length of a revolution (total FFT per period). In 
terms of geomagnetic geometry, the FFTs establish the duration for which the 
trajectory lies outside the trapping domain of the corresponding particle 
species, evaluated at the given energies. Or conversely, they are a measure 
of the degree to which the trajectory is exposed to the charged particle 
trapping domains. 
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The number of consecutive flux-free orbits of circular trajectories is 
primarily a function of altitude and inclination and to a lesser degree a 
function of particle energy. Of the investigated flight paths, and for the 
selected energies (electrons: E > .5 MeV, protons: FE > 5. MeV), none shows 
any completely flux-free revolutions per day. The total FFT, which derives 
completely from partially exposed revolutions (see “Exposure Analysis," Tables 
#31 to 36), in percent of total mission duration, can be summarized as 
follows: 


Altitude Protons Electrons 
(km) (E > 5. Mev) (E > .5 Mev) 
1667 22% 2% 

2593 24% 3% 
3889 30% 4% 
5186 35% 6% 
6389 39% 7% 
10371 51% 11% 


Higher energies will yield longer FFTs because the more eneryetic particles 
occupy a smaller volume of space. 


D. Dose and Shielding 


The calculated doses display features characteristic of the terrestrial 
radiation environment: at medium-to-high shield thicknesses, small 
contributions from relatively benign and low intensity electron spectra 
combined with major contributions from comparatively hard and intense proton 
spectra; at thin shield thicknesses the electrons predominate. 


Depending on altitude, the proton doses prevail for shield thicknesses from 
greater than 16 to greater than 200 mils of aluminum. At high altitudes (e.g. 
10371 km orbit) the opposite is true: protons are the major dose contributors 
only for shield thicknesses of less than 5 mils of aluminum. 


Significant, however, is the fact that for aluminum the proton dose is only a 
weak function of shield thickness, as it shows very little attenuation over 
the evaluated depth range. Thus, in order to get an appreciable reduction in 
the dose, say by one order of magnitude, a 2N-fold increase in shield 
thickness is necessary. The same is true for the bremsstrahlung dose. 
However, at low altitudes, in comparison to the proton contributions, the 
bremsstrahlung dose is so small (about 2 orders of magnitude lower) that it 
may be disregarded. 


I. Decay and Degradation 


The total doses obtained for each of the six investigated trajectories for the 
five year mission duration are substantial. In terms of electronics decay or 

materials degradation, the doses to be experienced on this mission inside the 

satellite, that is, behind a skin of about 80 mils of equivalent aluminum, are 
severe even for insensitive components or equipment: 


Altitudes (km) Spherical Dose 


1667 ~294K rads 
2593 ~1093K rads 
3889 ~1558K rads 
5186 ~1266K rads 
6389 ~1099K rads 
10371 ~1682K rads 


II. Contamination and Interference 


It should he remembered, that the direct or indirect effects of the radiation 
environment may also be a nuisance in terms of instrument interference or 
measurement contamination. If such is the case, some remedies may be 
available. 


III. Possible Improvements 


In the event that the magnitude of total dose or degree of radiation 
penetration behind the skin of the satellite is of importance to the mission, 
four possibilities exist to reduce the radiation effects on instrument and 
components: 


a) build or design an instrument less sensitive to radiation and 
construct the on-board and/or on-ground data processing software 
to remove or suppress radiation-induced noise 

b) change the orbit by any combination of the elements eccentricity, 
altitude, and inclination so as to achieve a more benign environment 

c) change the mission epoch: solar max for reduced proton intensities, 
solar min for reduced electron intensities 

d) provide increased shielding either by geometry or by weight or by a 
combination of both; 

by geometry: perform a 3-D analysis (solid angle sectoring) and 
rearrange other equipment on board the satellite in order to provide 
maximum protection to sensitive part over greatest possible fraction 
of solid angle.* 

by weight: place additional shields around sensitive part as needed. 
Clearly options (a), (b), and (d) are most readily accesible. 


*A powerful computer package for complex radiation shielding and transport 
calculations is now operational at GSFC. It is capable of addressing such 
topics as: (a) material mixtures, cross sections for protons, electrons, 
heavy charge particles, and neutrons, including source spectra and response 
functions; (b) source geometry, detector geometry, surfaces, rays, bodies, 
regions, body intersections, body unions, simple meshes, design bodies, 
spacecraft rays, with diverse features such as combinatorial options, 
translate-rotate-replicate capabilities, etc.; (c) heavy charged particle 
applications-1D transport by numerical integration, small volume pulse height 
(soft errors), 3D ray trace sectoring, 3D adjoint Monte Carlo; (d) electron 
bremsstrahlung-1D transport by numerical integration and hy adjoint Monte 
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Carlo, small volume pulse height, 3D ray trace sectoring, 3D forward and 
adjoint Monte Carlo, energy deposition, charging distributions. 
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